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Abstract 

The interaction of water molecules with metal surfaces is typically weak and as a result van der 
Waals (vdW) forces can be expected to be of importance. Here we account for the systematic poor 
treatment of vdW forces in most popular density functional theory (DFT) exchange-correlation 
functionals by applying accurate non-local vdW density functionals. We have computed the ad- 
sorption of a variety of exemplar systems including water monomer adsorption on Al(lll), Cu(lll), 
Cu(llO), Ru(OOOl), Rh(lll), Pd(lll), Ag(lll), Pt(lll), and unreconstructed Au(lll), and smah 
clusters (up to 6 waters) on Cu(llO). We show that non-local correlations contribute substan- 
tially to the water-metal bond in all systems, whilst water- water bonding is much less affected by 
non-local correlations. Interestingly non-local correlations contribute more to the adsorption of 
water on the reactive transition metal substrates than they do on the noble metals. The relative 
stability, adsorption sites, and adsorption geometries of competing water adstructures rarely differ 
when comparing results obtained with semi-local functionals and the non-local vdW density func- 
tionals, which explains the previous success of semi-local functionals in characterizing adsorbed 
water structures on a number of metal surfaces. 
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I. INTRODUCTION 



An accurate atomistic description of the water-solid interface is crucial for understanding 
many natural and technological processes such as atmospheric ice formation or fuel cell 
reactions. It is not surprising, therefore, that there has been an enormous amount of research 
devoted to this task in the past few decades^""-. In recent years, scanning probe techniques, 
in particular scanning tunneling microscopy (STM), have contributed significantly to the 
field by providing detailed insight into the structure and dynamics of water adstructures on 
the nanoscale. Generally such studies are, however, limited to well-defined, single crystal 
metal surfaces at low temperature and under ultra high vacuum (UHV) conditions''^. As 
a result the water-metal interface has become the workhorse system for understanding the 
basic chemistry and physics of how water interacts with solid substrates in general. 

Computer simulation techniques, in particular density functional theory (DFT), have also 
played a central role in understanding the structure of water on metal surfaces^i^. DFT has, 
for example, been essential in unravelling the structure of water overlayers on Pd(lll)^, 
Cu(llO)^, and Pt(lll)^. Notwithstanding the clear value of DFT in this area, the accuracy 
of the generalized gradient approximation (GGA) functionals, which are generally used in 
such studies, remains an important open question. This is mainly because the GGAs used 
fail to describe non-local van der Waals (vdW) dispersion forces, forces which are expected 
to be relevant in water adsorption and weak adsorption systems in general^'^'^'^. This casts a 
shadow over all standard density functional studies of water on metals and hampers progress 
in the field as it puts accurate predictions about interesting phenomena such as rates of 
heterogeneous ice nucleation beyond the reach of ab initio methods at present. 

Recent years have seen a number of exciting developments with DFT based schemes for 
dealing with vdW dispersion forces (see e.g. refs^^i"— and for a recent review^'). With this 
has come a flurry of interest in understanding the role of vdW forces in water adsorption 
systems. In particular, recent dispersion-corrected DFT studies have been reported for wa- 
ter adsorption on Rh(lll)25, Cu(110)26, Pd(lll)^>28^ Ag(lll)^, Au(lll)^'^, and Pt(lll)2a 
and by us on Cu(llO) and Ru(OOOl)— . These studies have indeed indicated that vdW dis- 
persion forces should be accounted for to properly describe the interaction of water molecules 
with the metal substrate underneath. In particular, it was shown that vdW dispersion forces 
explain a long standing discrepancy with respect to the relative thermodynamic stability of 
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bulk ice and wetting layers on Cu and Ru surfaces^. However, despite this work, important 
issues remain poorly understood, such as the relative importance of vdW forces on different 
metal surfaces. Another open question is the role of vdW in determining the most stable 
adsorption site and adsorption structures for hydrogen bonded water clusters on metals. 

Here we tackle these and related questions by exploring in detail the role of vdW disper- 
sion forces in the adsorption of water monomers and clusters on a number of metal surfaces. 
The particular approach we use to treat vdW is the non-local van der Waals density func- 
tional (vdW-DF) of Langreth and Lundqvist and co-workers^^ and some of its offspringi^. 
These functionals from the vdW-DF family have shown great potential when applied to 
a number of systems where dispersion forces are important^i22'2ii22ri^. From the current 
study we find that non-local correlations contribute substantially to the adsorption of water 
monomers and clusters on metal surfaces in general. The contribution to the binding from 
non-local correlations varies from substrate to substrate and is actually greater on the more 
reactive transition metal surfaces than on the noble metals. Although vdW enhances the 
bonding with the substrate, the water-water interaction energies within the clusters remains 
largely unaffected by the inclusion of vdW dispersion forces. In addition the geometries of 
the water monomers and clusters on the various surfaces are not affected to any great extent 
by the inclusion of vdW dispersion forces, which explains the previous success'^'^ of GGA 
functionals in predicting adsorption structures of water on metals. 

II. METHODOLOGY AND COMPUTATIONAL DETAILS 

Density functional theory and supercell periodic models were used within the VASP 
5.2 cod&-~'^'^. Total energies and electron densities were computed with various different 
exchange-correlation functionals: (i) the semi- local Perdew-Burke-Ernzerhof (PBE)^^; (ii) 
the non-local vdW-DF of Dion et al}'^; and (iii) a couple of modified versions of the original 
vdW-DF, where the original GGA exchange functional has been replaced by an optimized 
PBE (optPBE) or optimized Becke88 (optB88) in order to improve the accuracy of the vdW- 
DF schemei^. All vdW-DF calculations were carried out self-consistently within VASP as 
implemented by Klimes et al?'^ using the approach of Roman-Perez and Soler^^. In vdW-DF 
the exchange-correlation energy is calculated by adding three different terms: a GGA ex- 
change energy, the local correlation energy obtained within the local density approximation 
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(LDA) , and a non-local correlation energy based on electron densities interacting via a model 
response function. In the original vdW-DF the GGA exchange is obtained with the revPBE 
functional^. However in ref. p^] alternative exchange functionals to revPBE were developed 
which can improve significantly the accuracy of the vdW-DF method. In all calculations the 
core electrons were replaced by projector augmented wave (PAW) potentials'^ whilst the 
wavefunctions of the valence electrons were expanded in plane-waves with a cut-off energy 
of 600 eV. PBE-based PAW potentials were used for all calculations. A Monkhorst- Pack- 
grid with 12x12x1 k-point sampling per (1x1) unit cell was used. 

In the case of water monomer adsorption, the close-packed (111) and (0001) surfaces 
were modeled by p{2x2) unit cells, containing 6 atomic layers separated by at least 14 A 
of vacuum (23 A when computing binding curves). In the case of adsorbed water dimers 
on Cu(llO) we considered a p(3x5) unit cell containing 4 atomic layers, whilst in the cases 
of adsorbed trimers, tetramers, pentamers, and hexamers a p(4x6) unit cell was employed. 
The metal atoms in the 3 (2 in the case of adsorbed water cluster models) bottom layers were 
fixed at their bulk-truncated positions during structure optimization procedures. We note 
that PBE lattice constants were used for all functionals. This is likely to infiuence very little 
the adsorption energies obtained with vdW functionals, since, for example, lattice constants 
are typically within 0.010 of PBE for optB88-vdW^. In all cases a dipole correction along 
the direction perpendicular to the metal surface was applied^i^ and geometry optimizations 
were performed with a residual force threshold of 0.015 eV/A. 

Adsorption energies per water molecule were computed as follows: 

^[HsO/M] - E[M] - n£;[H20] 
n 

where i?[H20/M] is the total energy of the adsorbed n H2O molecules, E\M\ is the total 
energy of the relaxed bare metal slab and i?[H20] is the total energy of an isolated gas 
phase H2O molecule. Within this definition a negative adsorption energy corresponds to a 
favorable (exothermic) adsorption process. 

In order to analyze the role of van der Waals dispersion forces on adsorption, we have 
decomposed i^ads in to different energy contributions^. We define the water-water contri- 
bution, E^J, to Eads as: 
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where i?*°*[nH20] is the total energy of the water structure in the absence of the sub- 
strate, but with all atoms fixed in the precise geometries they adopt in the adsorption struc- 
ture. Notwithstanding the fact that any energy decomposition scheme is to some extent 
arbitrary^, the estimate of the water-metal bonding, -E^™, is simply taken as the difference 
between E'ads and E^^. We have also examined the non-local correlation part, E^^'^, of the 
exchange-correlation energy which is obtained directly from the calculation. We note that 
^nic inciij^fjgg a part which acts as a semi-local correction to LDA correlation, i.e. in a sim- 
ilar spirit as PBE correlation works. Thus the "non-local" interaction means interaction 
beyond LDA correlation. The corresponding attraction due to the non-local correlation in 
the adsorption energy per water molecule can be calculated as: 



an analogous expression to Eq. [H but where total energies are substituted by their non-local 
correlation part. 

III. RESULTS AND DISCUSSION 
A. Adsorbed water monomers 

First we consider the adsorption of an isolated water monomer on a series of metal 
surfaces. This is the simplest water adsorption system possible to investigate the role of 
vdW interactions with the metal, since no H-bonds between water molecules are present. 
In the first place we explored the effect of vdW interactions in determining the most stable 
adsorption site on a selection of close-packed metal surfaces. Considering the most stable 
adsorption site, we then extended the study to a larger number of substrates to investigate 
in a systematic manner the dependence of vdW interactions with the nature of the metal. 

Previous DFT studies with semi-local functionals have predicted that water adsorbs on 
atop sites on close-packed surfaces^"—. Since recent studies have shown that vdW forces can 
alter adsorption structurea^"^"— , it is important to establish if vdW changes the preferred 
site for water adsorption on metals. To this end we computed the optB88-vdW (and for 
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reference PBE) adsorption energies of a water monomer on Ag(lll), Au(lll), and Ru(OOOl). 
More than 30 different adsorption sites and molecule orientations were explored, from which 
6 stable and representative structures are shown in Fig. [TJ We found that the most stable 
adsorption site with the optB88-vdW functional is the atop site, consistent with previous 
work and our own new PBE calculations reported in Table [B These results indicate that 
accounting for vdW interactions does not lead to a qualitative change in the mechanism 
governing the water-metal adsorption geometry of adsorbed monomers^. Nevertheless, vdW 
interactions increase the strength of the bond with the substrate. In particular, on the noble 
metals, orientations of the water molecule with the oxygen atom away from the surface (S5 
and S6) have almost no binding to the surface at the PBE level, whereas optB88-vdW 
adsorption energies are relatively large (Table H]). This will very likely be of relevance to 
water monomer diffusion, with the diffusing water molecule able to access many more stable 
configurations when dispersion forces are accounted for. 

In order to gain more insight into the role of vdW forces on the adsorption energy 
and their dependence on the specific functional chosen, we calculated the PBE, revPBE- 
vdW, optPBE-vdW, and optB88-vdW binding curves of a water monomer on Ag(lll) and 
Ru(OOOl) as a function of the distance between oxygen and the metal atom underneath (Fig. 
[2]). Each point in the graph is obtained by keeping the z coordinate of the water oxygen 
atom fixed while relaxing the rest of the atoms in the system except the 3 bottom layers of 
the metal slab. Of the functionals considered, PBE predicts the weakest interaction at long- 
range. Although revPBE-vdW recovers the long-range attraction, the interaction strength 
around the equilibrium distance is only slightly larger than that obtained from PBE on 
Ag and on Ru the revPBE-vdW binding minimum is actually shallower than that obtained 
with PBE. At long-range, the energies of the optPBE-vdW and optB88-vdW functionals 
are similar to revPBE-vdW, but on approaching the surface become significantly more at- 
tractive. Computing the non-local correlation contribution to the total adsorption energy 
reveals that this is indeed the leading attractive term between the water molecule and the 
surface as observed previously by Hamada and co-workers^^. 

Considering the most stable adsorption site and adsorbed water orientation identified 
previously (SI), we examine now the adsorption energy of water monomers on a large range 
of metal substrates: Al(lll), Cu(lll), Cu(llO), Ru(OOOl), Rh(lll), Pd(lll), Ag(lll), 
Pt(lll), and unreconstructed Au(lll). Also included in this systematic study are results 
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from the original revPBE-vdW functional for comparative purposes with respect to optB88- 
vdW. The computed adsorption energies and optimized distances are summarized in Table 
Hn A key observation is that revPBE-vdW does not always enhance adsorption energies 
with respect to PEE. Indeed the revPBE-vdW adsorption energies can either be larger or 
smaller than PBE. This is consistent with previous studies for various adsorption systems 
where the revPBE-vdW adsorption energies can either be similar to or slightly smaller than 
PBFM^'^^^''—. This behaviour is a direct consequence of the underlying overly repulsive 
revPBE exchange in the revPBE-vdW kernel, which also causes pure revPBE adsorption 
energies to be underestimated compared to PBE. Indeed water-metal distances (dw-m) — 
defined here as the distance between the O atom of a water molecule and the nearest metal 
atom on the surface — optimized with revPBE-vdW are substantially larger (up to 0.32 A) 
than those obtained from PBE. Another key observation is that optB88-vdW consistently 
provides larger adsorption energies than PBE and similar water-metal distances (within 
0.1 A) for all investigated systems. As observed before in the binding curves of Ag(lll) 
and Ru(OOOl), we see that non-local correlation (-E"^s) is again the principal attractive 
contribution to the total adsorption energy (Eads) as shown in Table HIl Specifically, 
is in most cases considerably larger than Eads at the equilibrium geometry independent of 
the metal and functionals considered. Although the contribution of vdW dispersion forces 
to the non-local correlation at equilibrium distances is probably not dominant, it clearly 
enhances the total adsorption energy as indicated by the fact that typical optB88-vdW 
adsorption energies are in general 110-190 meV/H20 larger than PBE values (Table HTl) . 
Interestingly, larger adsorption energy enhancements are found on more reactive metals 
than on noble ones. For example, when comparing Pd(lll) and Ag(lll), the difference 
between PBE and optB88-vdW adsorption energies is 180 meV/Il20 and 144 meV/H20, 
respectively. Similarly, the vdW enhancement is also larger when comparing the adsorption 
on Pt(lll) and Au(lll): 186 meV/H20 and 158 meV/H20, respectively. This behavior can 
be rationalized in terms of the equilibrium water-metal distance in each CclSG clS discussed 
recently by Liu et al. for benzene adsorption on metal surfaces^. As shown in Table [III the 
water molecule adsorbs closer to the surface on reactive metals such as Pd and Pt than it 
does on more noble substrates such as Ag and Au. The shorter water-metal distances lead 
to larger vdW dispersion forces and, therefore, larger total adsorption energies. 

Overall we have shown that optB88-vdW (and for the systems considered optPBE-vdW) 
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increases adsorption energies by accounting for vdW interactions without substantially al- 
tering adsorption sites or bonding distances with respect to PBE. In contrast, although 
revPBE-vdW is able to catch similarly vdW interactions at long distances, adsorption ener- 
gies are not increased beyond those obtained with PBE. Unfortunately, experimental data is 
not yet available against which we can benchmark our computed adsorption energies. This is 
mainly because unlike simpler adsorption systems such as CO or noble gas adsorption, water 
molecules readily diffuse and cluster on the surface, forming complex cluster and overlayer 
structures (e.g. the clusters discussed in section ITlIB|) . Clusters and overlayers obviously 
complicate the analysis of experimental adsorption enthalpy determinations, especially with 
approaches such as temperature programmed desorption^. Although experimental adsorp- 
tion energies for water monomers on well-defined metal surfaces are lacking, we note that in 
a recent study it was possible to perform the first single crystal adsorption micro calorime- 
try measurements of the energy of a well-dened water adstructure (a partially dissociated 
water adlayer stabilized by hydrogen bonding) For this one overlayer structure we were 
able to show that the inclusion of vdW forces using optB88-vdW was crucial for achieving 
quantitative agreement with experiment^; a similar conclusion has also been drawn for the 
adsorption of benzene on Au(lll) and Pt(lll)^. However, in general, accurate experi- 
mental measurements of adsorption energies of water on well-defined surfaces are severely 
lacking and urgently needed^i^^. 

B. Adsorbed water clusters 

Let us focus now on more complex water adsorption systems than simple adsorbed water 
monomers, in which we simultaneously consider the interplay of water-metal and water- 
water interactions. To this end we have considered a series of adsorbed water clusters with 
between 2 and 6 water molecules on Cu(llO). We have considered Cu(llO) because of a 
number of recent STM experiments for low coverages of water on this surface^'^'^"— and 
because the open (110) surface allows us to explore a rich variety of isomers for a given 
water cluster and, in so doing, analyze the role of vdW in determining the relative energies 
of various adsorbed water clusters. 

For each of the clusters examined we considered a range of different initial geometries 
which upon optimization lead to a number of stable or metastable water clusters, as shown in 
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Fig. |3l In Table UTT] we summarize the adsorption energies of all these structures employing 
various functionals. For the adsorbed trimers and tetramers our results are in general agree- 
ment with a recent study by Kumagai et al?^^ who by combining STM and DFT identified 
chain trimers and cyclic tetramers as the most stable species. 

Consistent with the results of the water monomers we find for the clusters that including 
vdW interactions through optB88-vdW leads to a substantial increase (>120 meV/H20) 
of the total adsorption energy when compared with PBE. The optPBE-vdW adsorption 
energies are typically smaller than those from optB88-vdW by ^50 meV/H20 and revPBE- 
vdW values are in all cases less than PBE by at least 45 meV/H20. Although in most 
situations the structures obtained with PBE and optB88-vdW are similar, there are some 
occasions when accounting for van der Waals dispersion forces can alter the relative stabilities 
of the clusters considered. In particular, optB88-vdW tends to favor planar structures over 
buckled ones. For example, in the case of water tetramers the buckled cluster (Te-I) is 
preferred at the PBE or even revPBE-vdW level, over the planar tetramer (Te-II) by about 
10-15 meV/H20. However, the situation is reversed when considering optB88-vdW which 
predicts the planar tetramer (Te-II) to be 15 meV/Il20 more stable than the buckled one. 
It turns out that this result can explain the absence of corrugation in the STM images of 
Kumagai et al.—. The lack of buckling was attributed to tip-induced reorientation or fast 
dynamical fiuctuations between the two isomers^^. However, here we see that there is no 
need to invoke such arguments as the most stable structure when vdW is accounted for is 
planar. Similarly, the preference for planar adstructures when switching on vdW dispersion 
through optB88-vdW is again observed in the case of the P-II and P-III pentamers, and the 
H-II and H-III hexamers (Table [ml). 

Another important observation is that the clusters become more stable as the clusters 
increase in size (Fig. Hj). This is consistent with what has been observed on other surfaces^"^ 
and gas phase clusters21i2^ and is due to the formation of H-bonds between water molecules 
and to a well-known cooperative effect, i.e., the increase of the number of H-bonds within 
the clusters results in stronger H-bonds. Interestingly, the nature of this cooperative effect 
shows a very small dependence with respect to the particular functional considered: both 
PBE and the vdW-DFs present similar trends, only the magnitude of the total adsorption 
energy is shifted (Fig. H]). We see, therefore, that vdW forces seem to play a minor role in 
describing cooperative water-water effects in these systems. 
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In order to gain deeper insight into these adsorption systems, we decomposed the total 
adsorption energy into water-water and water-metal bonding contributions as given by Eq. 
[2J As shown in Table llllt the water- water bonding contribution to the adsorption energy 
obtained from the various functionals is very similar: the PBE, optPBE-vdW, and optB88- 
vdW results are all within 13 meV/H20 for all adsorption systems (the revPBE-vdW water- 
water interactions are systematically smaller in all cases). It is clear, therefore, that the 
main effect of vdW forces on the total interaction comes from the water-metal bonding. 
Indeed the approximately 120 meV/Il20 increase in the adsorption energy comes almost 
exclusively from the water-metal bonding. This is consistent with our earlier work on ice- 
like films on Ru and Cu^i. Unlike the GGA functionals where the contribution to the 
binding comes from density overlap, the non-local functionals give binding even between 
non-overlapping electron densities^. Therefore, when considering non-local functionals not 
only a small overlap region gives contribution to the adsorption energy, but a larger part of 
the surface too. This leads to a stronger water-metal interaction for non-local functionals 
than GGA functionals. Moreover, this also explains that the shift in the adsorption energy 
caused by the vdW density functional with respect to PBE is almost constant between 
different clusters. In general, water-metal bonding is greater than water-water bonding 
and is responsible for the relative stability of the isomers and preference for more planar 
rather than buckled isomers, independent of the considered functional. Upon comparing the 
strength of the water-water interaction for the various clusters considered, we find that the 
water- water bonding is strongest in the water tetramers (Te-I and Te-II). Larger clusters 
have larger total adsorption energies, but their water-water bonding is substantially reduced 
with respect to the optimal H-bond configuration offered by a tetragonal arrangement on 
the Cu(llO) surface. This becomes especially evident in the case of pentamers, where the 
most stable isomer, P-I, shows actually a relatively small optB88-vdW water-water bonding 
(-169 meV/H20) when compared for example to isomer P-V (-260 meV/Il20). 

Regarding the effect of vdW forces on the adsorption geometry of the clusters, we summa- 
rize in Table IIVI and Fig. |5] the averaged water-metal and water-water distances computed 
with different functionals. Essentially PBE geometries do not differ significantly with re- 
spect to the optB88-vdW or optPBE-vdW functionals for both water-metal and water-water 
distances. All distances are within 0.05 A of PBE for optB88-vdW and within 0.09 A of 
PBE for optPBE-vdW. On the contrary, revPBE-vdW predicts larger water-metal (<0.15 
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A) and water-water (<0.14 A) distances. It is interesting, therefore, that despite a noticeable 
enhancement in water-metal bonding when vdW is accounted for, the water-metal distances 
remain rather similar. 

IV. CONCLUSION 

We have investigated the role of vdW dispersion forces in water-metal bonding by consid- 
ering different non-local vdW-DFs and a range of water monomers and small water clusters 
on a series of metal surfaces. Analysis of our results reveals an enhancement of adsorption 
energies (typically >110 meV/H20) due to vdW interactions with respect to the widely 
used PBE functional. The increase in the adsorption energy comes almost exclusively from 
an increased water-metal interaction, with the water-water interaction within the adsorbed 
clusters being essentially unaffected by the inclusion of vdW dispersion forces. As a conse- 
quence, we observe that in general the explicit consideration of vdW dispersion forces does 
not alter the relative stabilities of structures predicted by PBE. In addition, despite increases 
in the total adsorption energies, the adsorption sites remain unchanged and adsorption ge- 
ometries (water-metal bond and water H-bond lengths) are very similar when comparing 
PBE and optB88-vdW in general. On the few occasions when the PBE and optB88-vdW 
structures differ, the optB88-vdW structures are flatter with any high-lying water molecules 
brought closer to the surface, which for the particular case of the water tetramer on Cu(llO) 
produces an adsorption structure in better agreement with experiment. This is a direct 
consequence of the fact that non-local correlation is enhanced by shorter water-metal dis- 
tances. It is interesting to note that the general similarity between the structures obtained 
with and without vdW dispersion forces is in contrast to what has been found for gas phase 
water clusters. In particular for the water hexamers — where this issue has been considered 
in greatest detail — the relative energies of the relevant low energy isomer structures differ 
completely depending on whether vdW is accounted for or not''^-. 
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FIG. 1: Top and side views of water monomers adsorbed at different sites on a close-packed metal 
surface. Small black, red, and grey spheres stand for H, O, and metal atoms, respectively. 
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FIG. 2: Adsorption energy of H2O monomer on Ag(lll) and Ru(OOOl) as a function of the vertical 
water- metal distance (dw-m), defined as the distance between the O atom of a water molecule and 
the nearest metal atom on the surface. Four different functionals are considered and the optB88- 
vdW non-local correlation (nlc) contribution to the total energy (Eq. [3]) is also shown (open circles 
and dashed line). The lines are merely a guide to the eye. 
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FIG. 3: Top view of water monomer (M), dimer (D), trimers (Tr-I and Tr-II), tetramers (Te-I and 
Te-II), pentamers (P-I, P-II, P-III, P-IV, P-V, and P-VI), and hexamers (H-I, H-II, H-III, and 
H-IV) adsorbed on Cu(llO). Small black, red, dark grey, and light grey spheres stand for H, O, 
and Cu in the first and second layer, respectively. Light red spheres indicate water molecules that 
are relatively far away from the metal surface (typically >3 A). 
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FIG. 4: Adsorption energies of different sized water clusters adsorbed on Cu(llO) using PBE, 
revPBE-vdW, optPBE-vdW, and optB88-vdW. Dashed lines connect the most stable isomers for 
a given number of H2O molecules with each functional. 
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FIG. 5: Averaged nearest neighbor water-metal, (top), and water-water, dw-w (bottom), 

distances for the water clusters depicted in Fig. [3] using PBE, revPBE-vdW, optPBE-vdW, and 
optB88-vdW. The lines connecting the points are there to guide the eye. 
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TABLE I: Adsorption energies (in meV/H20) of an isolated water monomer on three different 
metal surfaces with PBE and the optB88-vdW functionals. Six different adsorption sites and 
adsorbed water geometries are shown (see Fig. [T]). SI is the most stable adsorption structure on 
all surfaces with all functionals. Negative adsorption energies correspond to favorable (exothermic) 
adsorption. 





SI 


S2 


S3 


S4 


S5 


S6 


PBE 


Ag(lll) 


-137 


-134 
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-1 
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25 


Au(lll) 


-123 


-122 


-68 


30 
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62 


Ru(OOOl) 


-360 


-356 


-121 


-130 


-17 


-17 


optB88-vdW 


Ag(lll) 


-281 


-280 


-195 


-151 


-117 


-103 


Au(lll) 


-281 


-281 


-213 


-120 


-122 


-76 


Ru(OOOl) 


-541 


-538 


-271 


-334 


-131 


-125 
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TABLE II: Adsorption energies (-Eads in nieV/H20) and the optimized distance between the O atom 
of a water molecule and the nearest metal atom on the surface (dw-m in ^) for water monomers 
at the equilibrium adsorption site (SI in Fig. [T]) on all the metal surfaces investigated. Results for 



the PBhv, 


revPBh-vdW, 


and optB88-vdW functionals are reported. 






Metal 


PEE 


revPBE-vdW 


oDtPBE-vdW 


optB88-vdW 




-E'ads 


-E'ads 


-^ads 


-C/ads -C/ads 


-E-ads 


-^ads 


Al(lll) 


-192 


-165 


-380 




-302 


-476 


Cu(lll) 
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-201 


-360 


_9Q9 —ZL'^S 


-319 


-456 


Ru(OOOl) 
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-328 


-490 


-488 -530 


-541 


-563 


Rh(lll) 


-306 


-293 


-466 


-442 -527 


-497 
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Pd(lll) 


-243 


-263 


-423 


-380 -458 


-423 


-529 


Ag(lll) 


-137 


-192 


-260 


-267 -345 


-281 


-365 


Pt(lll) 


-217 


-241 


-391 


-358 -495 


-403 


-500 


Au(lll) 


-123 


-192 


-392 


-264 -398 


-281 


-433 


Cu(llO) 


-361 


-317 


-389 


-457 -486 


-503 


-474 




Al(lll) 


2.232 




2.503 


2.248 




2.195 


Cu(lll) 


2.507 




2.795 


2.495 




2.422 


Ru(OOOl) 


2.383 




2.505 


2.408 




2.361 


Rh(lll) 


2.359 




2.519 


2.389 




2.334 


Pd(lll) 


2.469 




2.714 


2.503 




2.419 


Ag(lll) 


2.788 




3.042 


2.746 




2.704 


Pt(lll) 


2.499 




2.815 


2.555 




2.498 


Au(lll) 


2.836 




3.084 


2.853 




2.750 


Cu(llO) 


2.185 




2.320 


2.202 




2.167 



22 



TABLE III: Computed adsorption energies (i?ads) of water monomer (M), dimer (D), trimers 
(TV), tetramers (Te), pentamers (P), and hexamers (H) adsorbed on Cu(llO) using PBE, rcvPBE- 
vdW, optPBE-vdW, and optB88-vdW DFs. Water-water (E^^) and water-metal {E^} bonding 
contributions are also given. All values are in meV/H20. 
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-361 


-317 
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-110 


-320 
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-113 


-307 


-564 
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Te-T 
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P-IV 
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-288 
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-339 
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-165 


-305 


-615 


-193 


-421 


-667 


-190 


-477 


H-III 


-523 


-193 


-330 


-463 


-166 


-297 


-613 


-195 


-418 


-668 


-192 


-476 


H-IV 


-514 


-156 


-358 


-456 


-141 


-315 


-607 


-163 


-444 


-661 


-155 


-506 
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TABLE IV: Averaged nearest neighbor water-metal ((i„_in) and H-bond ((i„_w) distances for 
different water clusters adsorbed on Cu(llO) using PBE, revPBE-vdW, optPBE-vdW, and optB88- 
vdW DFs. All values are in A. 
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z.ioo 


1 


77/1 


1/10 
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Z.Z 


1 VO/I 

i. ( y4 


P-I 


2.265 


1 
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i 
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1 80Q 
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1 
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1 
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2.208 


1 
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1 
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1 
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1 
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1 
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1 
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1 


801 
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1 
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1 
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1 
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1 
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1 
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1 
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2.102 
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1 
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1 
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1 
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1.805 
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2.143 


1 
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2.255 


1 


894 


2.164 


1 
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H-III 
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1 
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1 


926 
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1 


846 
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1 
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1 
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1 
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